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ABSTRACT 

A  general  method  for  the  analog  simulation  of  an  n   order  digital 
discrete  compensator  in  a  sampled-data  servo  system  is  presented,,  The 
method  employs  operational  amplifiers  such  as  are  found  in  ordinary 
laboratory  analog  computers 5  and  relay  sampling  switches  driven  by  a  low 
frequency  function  generator. 

The  method  is  discussed  in  relation  to  the  trend  to  the  use  of  dig- 
ital devices  in  control  systems  and  its  theoretical  development  in  the 
algebra  of  the  z-transformation  is  sketched .  The  validity  of  the  method 
was  proved  by  comparing  calculated  with  simulated  responses $  showing  the 
accuracy  of  the  simulated  system  to  be  no  less  than  that  of  the  computer 
itself e 

The  usefulness  of  the  simulated  discrete  compensator  is  shown  in  a 
general  investigation  of  step  responses  for  various  z-plane  root  loca- 
tions and  in  an  investigation  of  a  simple  adaptive  system. 


11 


ACKNOWLEDGEMENT 

The  author  wishes  to  express  his  appreciation  to  Lt.  Thomas  R. 
Ma  this,  U0  S0  Navy,  for  his  invaluable  counsel  and  assistance  with  this 
project.   In  addition,,  the  counsel  and  encouragement  of  Prof„  Richard 
Dorf  as  faculty  advisor  are  gratefully  appreciated. 

The  motivation  for  this  project  arose  during  a  wider  investigation 
into  a  digital-analog  control  system  conducted  by  U.  S.  Naval  Post- 
graduate School  officer  students  at  the  Dalmo-Victor  Company  of  Belmont,, 
California,  during  summer  I960.  The  cooperation  of  all  Dalmo=Victor 
personnel,  and  particularly  the  counsel  and  assistance  of  Mr.  Robert 
Mcintosh  and  Mre  James  Cox,  are  gratefully  appreciated. 


111 


iBLE  OF  Z   I 


Chai ' 


III 
IV 


Y 

IPI  5NDIX  I 

APPENDIX  II 

APPENDIX  III 

BIBLIOGRAPHY 


Title 
I   EDUCTION 
THEORETICAL  DEVELOPMENT 

2.1  System  Motivation 

2.2  Theory  of  Discrete  Compensator  Realisation 

2.3  Analysis  of  Delay  Circuit 

GENERAL  SIMULATION  OF  A  SECOND  ORDER  SYSTEM 

RESULTS  OF  INVESTIGATIONS 

'.I.I  Equipment 

h.2  Investigation  of  z-plane  Root  Locations 

Ii.3  Investigation  of  an  Adaptive  System 

CONCLUSIONS 

GRAPHICAL  EXPLANATION  OF  DELAY  CIRCUIT 

EQUIPMENT 

ANALOG  COMPUTER  CONFIGURATION  USED  IN 

INVESTIGATION 


IBM 
1 


10 
12 

l)i 
1)4 

V' 

■ 

39 
III 

1*3 
U5 


IV 


LIST  OF  ILLUSTRATIONS 

Figure                       Title  Page 

1  Typical  sampled-data  system  It 

2  Simplified  compensated  sampled-data  system  £ 

3  Realization  of  D.(z)  7 
U  Realization  of  discrete  compensator  9 

5  Block  diagram  of  sampled-data  system  8 

6  Analog  realization  of  time  delay  10 

7  General  simulation  of  second  order  system  13 

8  Uncompensated  system  z-plane  root  locus  l£ 

9  z-plane  root  locations  2\\ 

10  Step  responses  Roots  1  -  $  2% 

11  Step  responses  Roots  6-10  26 

12  Step  responses  Roots  11  -  1$  27 

13  Step  responses  Roots  16  -  20  28 
lit  Step  responses  Roots  21  -  25  29 
1$  Step  responses  Roots  26  -  30  30 

16  Step  responses  Roots  31  -  35  31 

17  Step  responses  Roots  36  -  UO  32 

18  Step  responses  Roots  Ul  -  U5  33 

19  Constant  peak  overshoot  loci  in  the  z-plane  3k 

20  Zero  migration  responses 

21  Pole  migration  responses  36 

22  Graphical  explanation  of  delay  circuit  U0 
23(a)  Sampler  and  hold  circuit  1*2 
23(b)  Relay  control  circuit  ll2 
2U  Analog  computer  configuration  used  in  investigation  UU 


BOLS  AND  ABBREVIATIONS 


Descri]  __^ 
a >   3!P  a10? * ' * *ai* ' *  * *am 


Coefficients  of  the  numerator  of  dis- 
crete compensator 


b,  b]  9   b-^QS. .  e5bi3.  ,.5b       Coefficients  of  the  denominator  of 

discrete  compensator 

C  H*  C(t)J  ,  C(s)?  C(z)       Output  signal 

c^(nT)s  c-j^  j_(n  -  l)TJ, ...  Output  time  sequence 

cps  Cycles  per  second 

D(z)s  D-,(z)s  Dp(z)  Transfer  function  of  discrete  compensator 

e^  Input  voltage 

e0  Output  voltage 

f  Frequency 

G(s)5  G(z)  Transfer  function  of  plant 

G^(s)s  Gjj(z)  Transfer  function  of  zero-order  hold 

circuit 

G^G(z)  Transfer  function  of  zero-order  hold 

circuit  and  plant  in  series 

k  An  integer 

K  Apparent  (root  locus)  gain 

K  Real  gain 

P 

K> 9   K, "  z-plane  gain 

m  An  integer 

Maximum  peak  overshoot 


n  An  integer 

p  z -plane  pole 

r(nT)  Input  time  sequence 

R  |"-  R(t)J  ,  R(s),  R(z)      Input  signal 

s  The  Laplace  variable 

vi 


Symbol 
T 
z 

■''. 

> 


Description 

Interval  between  samples 

z-transfcrm  variable  (  ~  e  ) 

z -plane  r,ero 

Second  order  plant  inverse  time 
constant 

Change  in  plant  inverse  time  constant 

Damping  factor 

Sampling  duration,,  sample  pulse  width s 
interrogation,  or  read  out  time 

Angular  frequency  (rad/second) 


f- 


I 


£ 


> 


b,   - 


Summer- 


Multiplication  by  a  constant  coefficient 


fr 


■ 


Inverter 


X* 


-        VT 


Integrator 


Sampling  sv.ri.tch 


Vll 


CHAPTER  I 


INTRODUCTION 


The  introduction  of  digital  computers  into  control  systems  has  led 
to  the  use  of  combinations  of  digital  and  analog  devices  within  the  same 
system,,  One  interesting  aspect  of  this  use  of  digital  devices  is  that 
of  the  analog  simulation  study.  How  does  one  go  about  simulating  a  dig- 
ital computer,  and  all  the  attendant  coding  and  decoding  devices  neces- 
sary for  its  use  in  the  control  loop,  with  the  analog  computer?  This  is 
a  question  that  the  preliminary  designer  faces  when  budgetary  or  equip- 
ment limitations  preclude  the  use  of  digital  components  directly. 

It  is  usually  necessary  to  treat  digital-analog  systems  within  the 
context  of  sampled-data  theory  (l,  2,  3).*  In  this  context 9   digital 
devices  are  not  too  difficult  to  simulate.  A  digital  computation <>  such 
as  multiplication  or  addition,  or  combinations  of  these  in  equations 9 
can  be  programmed  almost  as  easily  in  the  analog  computer  as  in  the  dig- 
ital. Digital  conversion  devices  such  as  shaft  encoders 9   diode  matrices, 
Shannon-Rack  decoders,  etc.,  can  be  simulated  with  unity  transfer  func- 
tions. These  simulations  are  simple  provided  the  digital  computation 
time  is  short  compared  with  the  interval  between  samplings. 

The  simulation  of  the  digital  discrete  compensator  by  analog  tech- 
niques is  not  so  simple,  however.  The  discrete  compensator  gets  its  name 
from  the  fact  that  it  operates  on  the  signal  only  at  each  sampling  in- 
stant. It  has  a  z-transform  which  is  a  ratio  of  n^   order  polynomials  in 
z.  Its  effect  is  similar  to  that  of  a  filter  compensation  in  the  continu= 

^-Numbers  in  parentheses,  e.g.,  (l,  2,  3),  refer  to  references  in  the 
Bibliography. 


0113  system.   The  difficulty  in  analog  simulation  lies  in  the  fact  that 
the  digital  realization  of  this  discrete  compensator  involves  storage  of 
signals  ,  that  is,  time  delay,  for  periods  of  from  one  to  n  sampling 
periods.  Such  storage,  while  elementary  in  the  digital  computer,  can  be 
achieved  only  indirectly  in  the  analog  computer. 

It  is  the  purpose  of  this  paper  to  present  a  simple  technique  for 
the  direct  analog  simulation  of  a   discrete  compensator  in  a  physical 
system  and  to  suggest  possible  extensions  of  the  technique.  The  paper 
is  divided  into  three  parts:  (l)  The  theoretical  and  practical  back- 
ground 5  wherein  the  detailed  motivation  for  the  technique  is  developed^ 
(2)  A  description  of  the  general  analog  simulation!  and,  (3)  A  discus- 
sion of  the  results  of  two  investigations  designed  to  test  the  validity 
and  usefulness  of  the  technique. 

The  directness  and  simplicity  of  the  technique  are  its  major  ad- 
vantages o  It  can  be  used  with  the  simplest  of  computers  for  a  wide 
variety  of  laboratory  or  classroom  problems.  Changes  in  sampling  rate^ 
compensator  coefficients,  or  plant  constants  are  effected  by  simply 
changing  potentiometer  settings. 

The  development  of  this  technique  was  a  direct  outgrowth  of  a 
project  involving  the  development  of  a  digital  comparator  for  use  in  a 
simple  servo  loop,  and  subsequent  digital  discrete  error  signal  compensa- 
tion of  a  two-phase  servo  motor.  Analog  simulation  of  the  discrete 
compensator  became  necessary  when  budget  limitations  permitted  the  dig- 
ital mechanization  of  the  comparator,  but  not  that  of  the  compensator. 


,PTER  II 

THEORST IG A  L  DEVELOPMENT 
2^,1  System  Motivation. 

A  typical  sampled-data  system  might  have  the  configuration  shown  in 
Fig.  1.  The  devices  shown  in  the  blocks  in  Fig.  1  are  typical  components 
that  might  be  encountered  in  practice.  Because  of  their  one  to  one  equi- 
valence  the  encoders  can  be  represented  by  sample  switches ,  the  decoder 
can  be  represented  by  a  zero-order  hold  circuit,  and  the  comparator  by  an 
error  summer.  The  block  diagram  then  reduces  to  that  of  Fig.  2(a) 0 

The  location  of  the  sampler  switches  can  be  moved  across  the  summer 
and  around  the  digital  compensator,  since  the  digital  compensator  only 
receives  and  operates  on  error  signals  at  the  sampling  instants.  The 
system  diagram  then  becomes  that  of  Fig.  2(b). 

This  is  the  classic  discrete  compensator  circuit  described  in  the 
literature.   In  arriving  at  it,  the  following  assumptions  are  made? 

(1)  All  encoders  are  interrogated  at  the  same  time,  and 
their  difference,  the  output  of  the  error  summer,  is  the  er- 
ror that  exists  at  t  =  nT  (where  T  =  the  sampling  or  inter- 
rogation period.) 

(2)  That  the  interrogation,  or  read  out,  time  (corre- 
sponding to  sample  pulse  width  "T  )  is  short  compared  to  T. 
This  permits  analysis  on  an  impulse  basis.   It  might  be 
pointed  out  that  this  is  generally  valid,  since  interroga- 
tion time  will  usually  be  of  the  order  of  a  few  microseconds, 
whereas  T  will  most  likely  be  longer  than  a  millisecond. 
Interrogation  time  might  be  longer  if  the  encoder's  output 

is  serial  binary,  but  this  is  treated  in  the  next  assumption. 
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FIGURE  2.     SIMPLIFIED  COMPENSATED  SAMPLED-DATA  SYSTEM 


(3)  '    bation  time  in  the  comparator  and  discrete  com- 
t  compared  with  T.  This  coi  ] 

if  less  recent  digital  logic  modules  are  used,  but  for  most 

modern  equipment  it  is  not. 

It  would  be  possible,  of  course,  to  introduce  the  compensation 
between  the  hold  circuit  and  the  plant.  This  would  be  very  close  to 
continuous  compensation.  However,  the  output  of  the  hold  circuit  is 
quantized,  and  if  the  assumption  is  made  that  the  digital  computer  must 
be  time  shared  (often  one  of  the  justifications  for  a  sampled-data  sys- 
tem), the  analysis  of  the  continuous  compensation  gets  quite  difficult 
and  its  generality  is  reduced.  In  addition,  it  is  often  true  that  the 
plant  is  airborne  and  the  computer  on  the  ground  with  a  digital  com- 
munication link  between  them  (another  justification  for  the  sampled-data 
system).   In  this  case  it  might  be  necessary  to  eliminate  all  excessive 
airborne  components,  thus  the  desirability  of  having  the  compensation  on 
the  ground,  in  the  computer. 
2.2  Theory  of  Discrete  Compensator  Realization. 

In  general,  the  digital  discrete  compensator  can  be  represented  by? 

an  +  aiz""1  +  apz"2  +  ...  +  ana~n 

D(z)  =  — 1 l ■ — D— 

bQ  +  b^z"-1  +  b2Z_<:-  +  ...-+  bnz-n 

n 

akz"k 


k=0 


n 

v- — v 

b0  +  /    bkz"k 
k-1 

The  only  restriction  on  the  generality  of  this  statement  is  that 
bQ  s  0.*- 

*Tou,  J.  T.j  Digital  and  Sampled-Data  Control  Systems,  McGraw-Hill  Book 
Co. j  New  York,  19h9 ,   pp  U32-3. 
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Ragazzini  and  Franklin*  demonstrate  the  realizability  of  D(z)   in  the 
following  manner: 


Let 


D(z)     = 


—1    .  -2  -m 

an  +•  a,z       f  aQz       +   ...    -fa   z 
Old  m 


12  v 

1  +  bxz"     +  bgz"'"  +  ...   +  bkz~K 


and  let  D(z)  contain  the  two  factors 


C(z) 
R(z) 


Di(z)  r 


1  +  bxz  X  f  b2z~2  +  ...  +  bkz~k 


-1 


-2 


D  (z)  =  a0  f  a1z"  +  a2z~  +  ...  + 


a  z 

m 


-m 


The  block  diagram  for  D-^(z)  is  shown  in  Fig.   3.      (h) 
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>x 


C,LCn-2)Tj  C,[f«-0Tj 


S+C  I 
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■ 


,(nT) 


/.  o 


Q..O 


Figure  3. 

The  z-transform  of  the  output  sequence  is  C]_(z)  and  the  n   pulse 
in  the  time  domain  is  c^(nT).  Each  time  delay  is  T  seconds.  Tracing 
the  signal  through  the  block  diagram  we  see  that 

Cl(nT)  =  r(nT)  -  blCl[(n-l)T]  -  b2Cl  [(n-2)T]  -  ... 

r(nT)   -  Cl(nT)  t  b1c1[(n-l)l]  +  b^  [(n-2)Tj  +  ... 


*Ragazzini,  J.  R.,  and  Franklin,  G.  F. j  Sampled-Data  Control  Systems,, 
McGraw-Hill  Book  Co.,  New  York,  19U8,  pp  76-7. 


Talcing  the  z-transform  of  both  sides,  there  results 


-1 


-2 


R(z)   -   G1(z)(l  f  blZ'x  f  b2z"^  +  ...  +  bkz"k) 
which  yields  D-^(z). 

The  total  output  C(z)  is  now  obtained  by  operating  on  C-,(z)  by  Dp(z)o 
c-,(nT)  is  the  present  output  sample  and  it  should  be  operated  on  by  Sq. 
Similarly,  c-^  (n-l)T  should  be  operated  on  by  a-, .  This  operation  on 
successively  earlier  signals  is  the  significance  of  z~m  in  the  D„(z) 
terms.  If  these  operations  are  then  summed  we  have  the  desired  D(z). 
This  is  accomplished  by  the  arrangement  in  Fig.  h. 

This  is  one  of  many  configurations  possible  to  realize  discrete 
compensation.*  It  has  the  advantage  of  requiring  a  minimum  number  of 
delay  elements,  and  for  this  reason  was  used  in  this  investigation. 

Such  a  compensator  could  be  programmed  reasonably  into  a  digital 
computer,  and  because  of  the  possibility  of  storage  between  sampling 
periods,  a  general  purpose  computer  could  readily  be  time  shared  to 
handle  this  and  other  problems. 

The  design  engineer  would  now  desire  to  know  how  a  given  compensa- 
tion of  this  sort  would  affect  his  system.  Analytically,  this  could  be 

obtained  by  taking  the  modified  z-transform 

R(z)  D(z)  G^  (z,m) 


C(z,m) 


1  +  D(z)  GhG  (z,m) 


of  the  system 
Q  CO 


♦*§> 


Dc*) 


i 

He   I 


Figure  £. 


*Touj  pp  UU9-U^U. 
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where  GTG"(z)  and  G,  G(z,m)  are  the  z  and  modified-z  transforms ,  respec- 
tively, of  the  zero  order  hold  circuit  and  plant  in  series. 

If  D(z)  is  of  high  order,  this  analysis  becomes  quite  laborious. 
If  a  digital  computer  is  available  for  analysis,  it  might  still  be  in- 
convenient to  vary  the  program  for  experimental  changes  in  the  a,  and 
b.  of  D(z).  An  analog  computer  simulation  immediately  suggests  itself. 

The  chief  problem  involved  in  setting  up  an  analog  simulation  is, 
of  course,  realizing  D(z).  Several  methods  have  been  described  in  the 
literature.  Sklansky  converts  D(z)  into  an  equivalent  R-C  network  in 
the  s-doniain  (£).  This  method  has  the  disadvantage  of  an  extremely 
laborious  conversion  of  zrplane  coefficients  (the  a.  and  b . )  into  work- 
ing coefficients  in  the  s-plane. 

The  difficulty  in  reproducing  a  direct  analog  of  D(z)  lies  in  de- 
veloping a  storage  and  delay  circuit.  Bigelow  achieved  this  by  using 
operational  amplifiers  with  sequencing  switches,  and  later,  dekatron 
tubes  (6,  7).  Dye  accomplished  this  delay  by  magnetic  amplifiers  (8). 
These  methods  all  have  the  disadvantage  of  requiring  extensive  and 
sophisticated  equipment. 
2.3  Analysis  of  Delay  Circuit. 

The  delay  circuit  used  in  this  development  was  obtained  from  the 
following  circuit: 


%{h  sre   G 


hold 


(s)  = 


Gh(s) 


Figure  6. 
1  -  e-Ts 
s 


and  G(s)  = 


1 

"Ts 
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The  z -trans form  of  the  closed  loop  is: 


eQ  GhG  (z) 


(2) 


ei         1  J-  GhG  (2) 


now, 

GhG(z) 

- 

1-z-1 
T 

(2   -  1) 

1 
77" 


Tz 


T2      (z  _  1)2      z  .  1 


then'  !2  (z)   =     1/(2-1)         a   1      m   __1 

e±  1  +  1  /  (z  -  1)     (z  -  1)  4-  1      z 


or  o         e.  s  e.z 
?  0      1 

which  means  that  the  output  voltage  is  the  same  as  the  inpui  ,  but  delayed 
by  T  seconds „  This  is  true  provided  the  output  is  sampled^  for  the  above 
is  true  only  at  the  sampling  instants.  A  graphical  explanation  of  the 
above  phenomenum  is  contained  in  Appendix  I. 

This  delay  circuit  is  simple  to  set  up  on  any  analog  computer j 
requires  no  expensive  or  even  unusual  equipment;  and  is  adaptable $   since 
a  change  in  sampling  rate  requires  only  one  compensatory  potentiometer 
setting  change* 
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CHAPTER  III 

GENERAL  SIMULATION  OF  A  SECOND  ORDER  SYSTEM 

,  7  shows  an  analog  simulation  of  the  system  of  Fig.  2(b)  that 
is  valid  for  any  order  D(z)  where 

D(z)   = 


and         G(s) 


s  ( s  f  cA  ) 

with  sampling  frequency  =  l/T. 

This  is  an  exact  analog  of  the  digital  discrete  compensation  of 
Fig.  !u  There  is  no  explicit  storage  device  in  the  analog  simulation,, 
Since  analog  information  is  available  instantaneously ^  this  is  unneces- 
sary o  However,  because  this  analog  information  is  available  at  all 
times 3   tut  is  correct  only  at  t  =  nT  (the  sampling  instants )s  there 
exists  the  requirement  that  the  output  of  the  compensator  be  sampled $ 
and  that  all  samplers  be  in  synchronism.  This  latter  requirement  is 
easily  met  by  using  the  same  triggering  source  for  the  samplers,, 

In  practice  this  configuration  would  be  raoiified  considerably,, 
since s   in  general 9   the  plant  would  be  more  complex  than  second  order 
and  the  D(z)  coefficients  would  have  positive  or  negative  signs.  Sign 
changing  is  accomplished  routinely  by  inserting  sign  changing  amplifiers 
where  necessary. 
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FIGURE  7.  GENERAL  SIMULATION  OF  SECOND  ORDER  SYSTEM 
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CHAPTER  IV 

RESULTS  OF  INVESTIGATIONS 

hd  Equipment. 

The  analog  computer  used  was  a  Donner  Model  3000.  Tito  sampler  and 
hold  circuits 9   each  containing  a  relay  sampler  triggered  by  a  low  fre- 
quency function  generator  and  a  Philbrick  USA-3  operational  amplifier 
were  used,  A  more  detailed  description  of  these  is  contained  in 
Appendix  II, 
It  .2  Investigation  of  z -plane  Root  Locations. 

To  evaluate  the  validity  and  usefulness  of  the  technique 9   a  simple 
system  was  postulated  and  various  first  order  compensations  were  applied 
to  it.  This  permitted  a  relatively  simple  analytical  check  of  the  ac- 
curacy of  the  system, 

A  plant  with  the  Laplace  transform 
G(s)  =  — ~— 

s  (s  +  l) 

was  set  up.  A  sampling  period  (T)  of  one  second  was  chosen5  and  K  -  lo0 
was  the  initial  gain. 


1   e"Tj 
With    Gh(s)   =  =-^— 


,    s            (1  -  e   ) 
then        GhG(s)  -     - ~- 

s2(s  +-  1) 

—      0.368  (z  +  0,718) 
and        GhG(z)  "  (z  -  1)  (z~-  0.168) 

The  closed  loop  roots  for  this  system  with  no  compensation  (sampled  error) 
were  found  on  the  z-plane  root  locus  diagram  in  Fig,  8.   If  the  pole  at 
z  "  0.368  is  cancelled  by  a  zero  of  the  compensator 9   and  another  pole  is 
introduced  by  the  compensator,  by  changing  the  compensator  pole  and  the 

Hi 


system  gain,,  it  is  possible  to  produce  two  closed  loop  roots  anywhere  in 
the  z-plane.  This  was  done  for  most  of  the  stable  region  of  the  z-plane 
with  results  discussed  below. 

Table  I  and  Fig.  9  show  the  compensated  root  locations  for  the  var- 
ious compensations  investigated.  The  system  has,  for  a  first  order  com- 
pensator using  pole  cancellation,  the  open  loop  transfer  functions 

D(z)GhG(z)     (z  _  fe)       (z  _  Qi368)  (z  _  .±) 

The  roots  are  enumerated  and  identified  by  b  and  K'  in  Table  I  and 
located  in  the  z-plane  in  Fig.  9.  Note  that  each  complex  root  in  Fig.  9 
has  a  complementary  root  in  the  negative  imaginary  region,,  Step  re- 
sponses for  all  root  locations  are  contained  in  Figs.  10  -  18 „  The 
responses  for  root  locations  1,  2,  18,  and  20  were  compared  against 
analytically  computed  curves  with  results  that  corresponded  with  an 
error  of  %%   or  less.   It  is  felt  that  this  is,  in  general,  sufficient 
accuracy  for  the  design  engineer  doing  preliminary  work,  and  it  is  all 
that  can  be  expected  from  the  computer  without  chopper  stabilized  ampli- 
fiers and  temperature  stabilized  components.  In  any  event,  these  com- 
parisons with  the  relatively  laboriously  calculated  curves  are  sufficient 
justification  for  the  approach,  since  they  show  that  a  good,  quick  solu- 
tion of  the  digitally  compensated  plant  response  can  be  obtained  for  all 
time,  not  just  the  sampling  instants,  by  analog  simulation. 

Compensating  in  the  z-plane  doesn't  offer  as  many  guide  lines  as  does 
compensating  in  the  s-plane.  Among  other  things,  constant  V  lines,  which 
map  into  the  z-plane  as  logarithmic  spirals,  lose  their  meaning  there  due 
to  intersampling  effects.  Jury,  in  discussing  design  parameters,  takes 
as  one  such  parameter  the  maximum  peak  overshoot,  1L.,  and  draws  curves 
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for  this  in  the  z-plane.*  The  investigation  confirmed  Jury's  curves  and 
Ln  Lcated  that  Mq  might  be  the  best  simple  parameter  for  compensation „ 
Tnese  curves  are  extended  in  Fig.  19  for  most  of  the  z -plane. 

It  can  be  seen  in  Table  I  that  there  is  a  reasonable  correlation 
between  expected  transient  frequency  and  that  observed,,  although  this 
fact  is  of  questionable  usefulness  as  a  design  parameter.  Nothing  could 
be  deduced  from  the  response  data  that  served  to  link  root  location  and 
settling  time.  In  general  the  step  response  becomes  more  heavily  damped 
as  the  roots  move  toward  the  positive  real  axis  (0  -  z  -  1.0). 

Some  interesting  responses  are  produced  by  roots  on  the  real  axis. 
The  presence  of  a  negative  real  root  causes  ripple  even  where  tnere  is  a 
definite  positive  real  root  (root  location  #6).  This  is  true,  as  well, 
in  those  cases  where  the  residue  of  the  negative  real  root  vanishes,  e.g., 
the  case  of  cancellation  compensation  (//s  7,  8,  9).  This  ripple  can  even 
cause  overshoot  and  give  the  impression  of  underdamping  ( /9,  10). 

Root  #11  was  placed  at  the  stability  limit,  and  produced  the  expected 
erratic  phenomena  with  the  system  nonlinearities  taking  charge.  Among 
these  phenomena  was  the  marked  undershoot  following  the  first  overshoot. 
This  was  also  noted  for  the  complex  roots  located  in  the  vicinity  of 
z  -  -  1.0  (#s  3ii-36,  39-140  ♦  Roots  in  the  remote  stable  portions  of  the 
second  quadrant  near  the  unit  circle  display  stable  characteristics,  but 
have  far  greater  overshoots  and  undershoots  than  do  comparable  root  loca- 
tions for  continuous  systems. 

Roots  jhx^>   consist  of  a  stable  positive  real  root  and  an  unstable 
negative  real  root.  The  resulting  limit  cycle  demonstrates  this 


*  Jury,  E.  I.j  Samoled-Data  Control  Systems,  John  tfiley  &  Sons,  New  York, 
1958,  pp  132-U. 
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instability, 

resp  Le. 

kgA.JL    '  "-  ; ' ' '  an  A      e   'System, 

i  j  '■  e  defined  herein  to  be  a  s 

itself  to  produce     bter  response  as  certain  extern:]     lit: 
change  than  the  same  system  would  without  the  adaptive  feature.  In 
order  to  utilize  a  digital  computer  to  effect  changes 
adapi ive,  programming  in  advance  would  probably  hive  to  be  employed. 

an  example  of  this  use  of  the  digital  computer s   consider  the  f 
lowing  systems  Assume  that  it  is  known  that  the  time  constant  of  the 
plant  would  change  during  operation  due  to  environmental  changes s   such 
are  encountered  by  an  aircraft  or  missile  operating  under  varying  veloci- 
ties.,  densities,,  and  temperatures,  and  that  these  changes  occur  between 
well  defined  limits.  These  changes  can  be  measured  externally  to  the 
control  loop.  Assume  that  the  amount  of  change  of  the  time  constant  is 
a  definable  function  of  the  measured  parameters s   and  its  formulation  can 
be  programmed  into  a  computer.  Once  formulated,  the  time  constant  change 
is  then  applied  to  the  digital  compensator  to  effect  either  compensator 
pole  or  zero  shifts  or  both  to  compensate  for  the  system  pole  and  zero 
shift. 

To  investigate  a  system  of  this  type  it  '/as  assumed  that  the  sysr  ■ 

was  compensated  by  a  first  order  compensator,  D(z),  where 

U+  aio) 
D(z)  =  ______ 

(z  +  blc) 
arc  '    this  compensation  was  satisfactory. 


) 


1  ,  s  the     •  r  function 


s    (s    +         ) 


se  of  the  plant  time  constant.,    oc  3   is   varied  an  amount  £   A 
This  change  manifests  itself  in  the  z-transform  by: 

(1)  a  new  pole  position 

(2)  a  new  zero  position 

e_(a±  AoC)T  (T  4.  D  _  ! 

21  "  T+e-(XiZ-Vv)T^l 

(3)  a  new  z-plane  gain 

The  designer  has  at  his  disposal  a-i  and  b-,  to  re-compensate  for  the 
changed  conditions „  If  some  formula  of  the  form  either 

al  "  f^a105  b13>     ) 
or     b1  -  g(a103  b10,  .:  v  ) 

or  both s  can  be  determined  that  will  satisfactorily  compensate  the 
changed  system^  this  formulation  could  be  readily  programmed  into  the 
digital  computer.  This  can  be  called  pole  and  zero  migration,, 

In  the  investigation^  an  empirical  approach  to  obtain  such  formulas 
was  used0  The  system  with  initial  compensation  was  set  up  on  the  com- 
puter and  was  varied  an  amount  :Vi.  Py  trial  and  error 9  compensator 
pole  or  zero  migrations  were  introduced  and  the  step  response  observed „ 
The  criterion  for  success  was  the  minimum  deviation  of  the  step  response 
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from  the  or  pensated  response  in  overshoot  and  rise 

For  the  system  investigated: 

K  =  loOO  (K«  -     0,36°-) 
*     =  1.00 
a10  -  -0.368 
b10  =  -0.100 

This  was  the  compensated  system  with  root  location  '/3S  described  in  the 
previous  section.  Rise  time  was  satisfactory  with  Mp  =  18$. 

Zero  Migration  Qnly._ 

Satisfactory  response  was  demonstrated  for  -  0.6  -  /_        ~   0.2  when 

91  =  ^a10  +  Oo7   v  v  ^ 

Pole  Migration  Only. 

Satisfactory  response  was  demonstrated  for  -O.U  -      -  0.0  vriien 

bx  =  (b10  -  A  >c ) 

These  responses  are  displayed  in  Figs.  20  and  21. 

Although  a  system  could  be  envisioned  with  the  transfer  function 

G(S)   =  _5 _ 

s(s  4-  OL)  (1) 

a  more  likely  physical  second  order  system  (servomotor)  would  have  the 
transfer  function 

Kp  Ko 

G(s)  -       /        ,  -^     =     _™^-_ 

S  (1/oc  s  f  1)  s  (s  r     <-) 

where  K  is  the  real  gain  of  the  system.  K  in  (1)  above  is  now  a  func- 
tion of  the  plant  time  constant.  For  example ,  if  °i   were  decreased 5  the 
system  would  tend  to  become  more  oscillatory 3  but  simultaneously^  the 
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<-K_,  se  time        ".crease, 

nvesti  ga    \9   rise  bj 

;  for  r 
■  sator  alor     rase  this  compensator  c        0]         ~he 
error  signal  until  T  seconds  after  the  step  is  applied.  Rise  time  is 
largely  determined  by  the  slope  of  the  initial  rise  and  this  in  turn  is 
determined  by  the  system  gain,   2CKp.  If  T  were  very  much  shorter  than 
rise  time  this  kind  of  compensation  would  be  possible 9   but  then  the 
sampled-data  system  would  approach  being  a  continuous  system0 

However 9   if  the  gain  of  the  system  is  multiplied  by  the  new  time 
constant ^  1/ X-,  3   where  X     ■»  .X.  i  ZNOd  9   as  the  system  changes 9   then 
the  open  loop  plant,  transfer  becomes 

G(s)   =   (1/  d)  ^-P 

s(s  +   Xn) 

.7  ich  is  the  original  system  investigated  with  K  =  Kp  -  constant,,  This 
programming  of  the  gain  can  be  mechanized  easily 5  utilizing  the  same  in- 
formation source  from  which  Z\<X    was  obtained. 

It  should  te  noted  that  responses  obtained  for  pole  or  zero  migra- 
I  ions  are  satisfactory  for  a  limited  range  of   zi'VX.  More  elaborate 
re-compensation  by  formula  by  combinations  of  pole  and  zero  migration 
is  possible  but  was  not  investigated. 
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FIGURE  10.   STEP  RESPONSES  ROOTS  1-5 


FIGURE  11.  STEP  RESPONSES  R       -  10 
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FIGURE  12.     STEP  RESPONS 
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FIGURE  13.  STEP  RESPONSES  ROOTS  16  -  20 
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FIGURE  lh.     STEP  RESPONSES  ROOTS  21  -  2£ 
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FIGURE  15.     STEP  RESPONSES  ROOTS  26  -  30 
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FIGURE  16.     STEP  RESPONSES  FOOTS   31  -  3£ 
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FIGURE  17.     STEP  RESPONSES  DOOTS  36  -  1:0 
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FIGURE  20.   ZERO  MIGRATION  RESPONSES 
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valid  Lne  il        tion  of  a  disc        >n- 

sat  led-data  system  such  as  might  be  encountered  when  a  dig- 

ital computer  is  introduced  into  the  servo  loo  <,  Other  investigations 
of  interest  were  the  exploration  of  responses  for  various  root  locations 
in  the  z-plane,  and  a  brief  study  of  a  simple  adaptive  system  utilizing 
the  discrete  compensator.  The  conclusions  may  be  summarized  as  follows t 

(1)  It  is  completely  feasible  to  build  an  analog  simulation  of  a 
discrete  compensator  of  any  order  u4  li   ^    : nary  laboratory  com- 
ponents ,  the  order  of  the  compensator  being  subject  to  the  size  of  the 
computer  available „  The  limits  in  accuracy  are  commensurate  with  those 
of  the  computer ,  and,  in  any  case,  are  satisfactory  for  preliminary  de- 
sign work ,  averaging  less  than  $%   in  this  inv     iti<     If  1  e  digital 
computer  being  simulated  has  appreciable  computation  time,  this  can  be 
simulated  by  using  a  diode  bridge  sampling  switch  triggered  by  a  puis 

Layed  by  the  computation  time. 

(2)  The  simplicity  of  investigation  of  a  wide  variety  of  z~] 

roots  without  the  labor  of  modified  z-transform  calculations  suggests  an 
unlimited  usefulness  for  the  technique  in  practical  and  theoretical  in- 
vestigations into  optimum  compensation,  intersampling  ripple,  self- 
adaption,  etc. 

(3)  If  the  compensated  system  time  constant  changes  due  to  external 
causes,  and  this  change  is  known  to  be  within  certain  limits  and  calcu- 
able  from  externally  measured  parameters,  it  might  be  possible  to  effect 
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: 

by  its  zeroes.  If  thf  La  ted  were  a  ,  then  as 

the  motor's  time  constant  changed,  so  would  the  apparent  gain.  ] 
tern  coul.i  be  re-compensated  by  the  same  formula  mentioned  above  if  the 
gain  were  programmed  to  remain  constant  as  the  time  constant  changed. 

The  investigation  of  the  adaptive  system  was  quite  limited^  and  the 
possible  usefulness  of  this  method  suggests  that  further  investigation 
is  warranted. 
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PLANATJ      DELAY  CIRCUIT 

The  circuit  of  Fig.  6  can  be  considered  to  f miction  in  the  follow- 
ing manners  Suppose  there  is  a  variable  input,  signal  such  as  shown  in 
Fig.  22,  Since  at  zero  time  there  is  no  feedback  signal^  the  error  is 
the  same  as  the  input  and  the  integrator  starts  integrating  at  the  rate 
controlled  by  the  error  voltage  (supplied  by  the  sampler  at  t  -   0) 9   such 
that,,  assuming  a  linear  integrator 5  at  T  seconds  the  integrator  output  is 

T  x  (Error  Voltage)  x  l/T  =  Error  Voltage 
The  feedback  serves  to  progressively  eliminate  the  error  voltage.  This 
can  be  seen  at  time  t  =  T,  where  the  sampled  input  voltage  is  the  same 
as  it  was  at  t  =  0$   the  feedback  voltage  is  also  the  same  magnitude  but 
negative s  and  the  error  voltage  goes  to  zero.  The  integrator  holds  its 
voltage  with  no  error  signal,  and  commences  to  integrate  again  at  t  -  2T5 
when  it  receives  a  new  error  signal.  The  integrator  output  sampled  is 
thus  seen  to  be  the  image  of  the  input  signal  displaced  one  sampling 
period. 
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ir  Scientific  Co.  Model  3000  electronic  analog  computer  is 
a  semi -port able 5  ten  amplifier,  unstabilized,  relatively   i  .     Lve 
computer  which  is  generally  used  for  el     jry  classroom  and  laboratory 
work.  The  computer  has  reasonably  good  characteristics  for  its  si%e  and 
cost.  Each  operational  amplifier  is  a  stable,  high  gain  D.C.  circuit,, 
with  a  pentode  driving  a  cathode  follower.  Average  gain  is  better  than 
10,000,  Input  impedance  is  that  of  an  open  grid  pentode,  and  output 
impedance  is  less  than  one  ohm.  Long  term  drift  is  -  limv/hr.  Output 
range  is  =100  to  flOO  volts,  with  load  currents  up  to  £ma. 

The  Donner  amplifiers  were  not  satisfactory  for  use  in  th< 
and  hold  circuits,  however,  apparently  due  to  excessive  loadi 
this  purpose  the  chopper  stabilized  Philbrick  Model  USA  -  3  was  used. 
Plotters . 

All  responses  were  recorded  initially  on  the  F,  L„  Moseley  Co0 
P0GRAF  Model  23  X-Y  plotter  running  in  the  time  sweep  mode.  This 
>rder  permits  measurements  accurate  to  1  part  in  1,000. 

The  responses  contained  in  this  report  were  recorded  on  the  Clevite 
Go,  BRUSH  dual  channel  recorder  and  amplifier  Model  RD5621, 
Sa mpler  and  Zero-order  Hold  Circuit, 

Figs.  23(a)  and  (b)  contain  the  circuit  diagrams  of  the  sampler  and 
zero     r  hold  circuit  and  the  relay  control  circuit, 
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FIGURE  23(b).  RELAY  CONTROL  CIRCUIT 
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The  sampler  and  hold  circuit  is  described  in  Appendix  II,   In 

order  to  account  for  the  sign  of  a  and  b  3   the  sign  changing  ampli 
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fier  '/£  is  either  in  or  out  of  the  circuit  according  to  the  following 
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FIGURE  2U.  ANALOG  COMPUTER  CONFIGURATION  USED  IN  INVESTIGATION 
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